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Quantum disordered problems with a direction (imaginary 
vector-potential) are discussed and mapped onto a superma- 
trix (7-model. It is argued that the OD version of the a-model 
may describe a broad class of phenomena that can be called 
directed quantum chaos. It is demonstrated by explicit cal- 
culations that these problems are equivalent to problems of 
theory of random asymmetric or non-Hermitian matrices. A 
joint probability of complex eigenvalues is obtained. The frac- 
tion of states with real eigenvalues proves to be always finite 
for time reversal invariant systems. 
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New very interesting phenomena may occur in sys- 
tems with non-Hermitian quantum Hamiltonians. Al- 
though the Hamiltonian of a closed system in equilibrium 
must be Hermitian, non-Hermitian models can describe 
nonequihbrium processes ||l[ , open systems connected to 
reservoirs B , dynamics of neural networks H , and have 
many other applications. 

In a recent remarkable work Q Hatano and Nelson 
considered a model of particles described by a random 
Schrodinger equation with an imaginary vector-potential. 
This model arises as a result of mapping of flux lines in 
a (d + l)-dimensional superconductor to the world lines 
of d-dimensional bosons. Columnar defects introduced 
experimentally in order to pin the flux lines M lead to 
the random potential in the boson system, whereas the 
component of the magnetic field perpendicular to the de- 
fects results in the constant imaginary vector-potential 
ih [pi. The prediction made in Ref. B is that already 
a one-dimensional chain of the bosons has to undergo a 
localization-delocalization transition as a function of dis- 
order. This effect is new and very unusual for physics 
of disordered systems, which clearly shows that random 
models with the imaginary vector potentials deserve a 
further discussion. 

The Hamiltonian H of the simplest model of non- 
interacting particles with a disorder and the imaginary 
vector potential ih has the form 



H = ip+ihf I (2m) -|- U (r) 



(1) 



where p — — iV and U (r) is a random potential. 

Although the imaginary vector-potential i\\ appears 
quite naturally in the model of the flux lines, the exis- 
tence of such a quantity in the conventional quantum me- 
chanics seems to be rather exotic because it violates the 



hermiticity. However, one can imagine situations when 
this term may be finite. 

Of course, this is not possible for a closed system of 
the particles. But we know that taking into account 
degrees of freedom of an environment can lead to the 
non-hermiticity. One of the well known examples is a 
quantum dot connected to leads. Instead of considering 
the whole system one can eliminate the leads (see, e.g. 
||7|). This results in complex eigenenergies in the dot; the 
imaginary parts of the eigenenergies describe the proba- 
bility of an escape into the leads. 

At the same time, the non-Hermitian Hamiltonian, Eq. 
(|l|), differs from the Hamiltonians for open dots. It is 
real and contains the vector h that introduces a direc- 
tion. Average physical quantities have to depend on the 
direction of h but the time-reversal symmetry is not bro- 
ken. To understand better the physical meaning of the 
vector h it is instructive to write a lattice Hamiltonian 
Hl corresponding to H 
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where c"*" and c are creation and annihilation operators, 
and {g.vS are the unit lattice vectors. Considering the one 
particle Hamiltonian we do not need to specify statistics 
of the particles. (Eq. (g) was used in Ref. Q for numer- 
ical calculations). 

In Eq. (Q), the hopping probability along h is higher 
than in the opposite direction. In other words, the Hamil- 
tonian iJi describes a directed hopping in a random po- 
tential. This model can be considered as a quantum 
counterpart of a directed percolation model introduced 
by Obukhov S. Such models may describe properties 
of disordered systems without a center of inversion. The 
finite h in Eq. (g) can be due to, e.g. a possibility of 
a virtual tunneling into a dielectric in an electric field. 
This would not lead to an imaginary contribution to the 
chemical potential as in open quantum dots but might 
make hopping in different directions nonequivalent. It 
is interesting to note that Green functions used in Ref. 
B contained as the energy spectrum the combinations 
(p-|-ia) with a constant vector a, which corresponds to 
the continuum version of the Hamiltonian, Eq. (H^) . 

"Conventional" (non-directed) disordered systems ex- 
hibit a variety of different phenomena. Depending on 



the geometry of the sample and strength of disorder one 
can have conduction, locahzation or, e.g. quantum chaos. 
One can guess that the models described by Eq. (|]j|) also 
contain these effects although the effects may be peculiar. 
In analogy with the directed percolation it is reasonable 
to call the corresponding phenomena directed quantum 
chaos, directed localization, etc. 

The main results of Ref. B were derived for ID and 
2D samples from numerical computations, although some 
qualitative conclusions were made considering a one- 
impurity model. At the same time, very well developed 
analytical methods of study of conventional disordered 
systems exist and it is highly desirable to develop anal- 
ogous schemes for the directed disorder problems. This 
would help to understand the localization-delocalization 
transitions but also consider new phenomena like directed 
quantum chaos. 

In this Letter it is shown that problems of directed 
disorder described by Eqs. (|]j|) can be reduced to 
calculations within a non-linear supermatrix cr-niodel. 
This approach proved to be very useful for a broad va- 
riety of different problems (for a review see H). The 
cr-model derived below is applicable in any dimension 
and differs from previous ones by a new h-dependent 
term. Problems of quantum chaos correspond to the zero- 
dimensional version and are most simple for calculations. 
Some new results are obtained for this case classified be- 
low as directed quantum chaos. 

Recently a "ballistic" cr-model was derived by averag- 



ing over either rare impurities or energy 1 11 1. Apparently, 
proper h-tcrms can be written for this case, too. 

Although the imaginary vector potential h enters Eqs. 
(|l|,p|) almost in the same way as the physical vector- 
potential A, the presence of i changes considerably 
the derivation of the tr-model. A simple replacement 
(e/c) A ^ ih in the cr-model of Refs. |l9|,|l0| , would lead to 
the absence of the ground state. This is because eigenen- 
ergies of the Hamiltonians, Eq. (|l]) or Eq. (||), are not 
necessarily real and retarded Gf and G^ Green functions 
can have now poles everywhere in the complex plane of 
e. As a result, these functions cannot be written in the 
usual form of convergent Gaussian integrals, which was 
a necessary step in derivation of the cr-model. 

First of all one should choose a proper quantity to cal- 
culate. According to a discussion of Ref. Q an important 
information can be obtained studying a joint probability 
P (e, y) of real e'j, and imaginary e^. parts of eigenenergies 
or, in other words, the density of complex eigenenergies. 
This function is introduced as 
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where V is the volume of the system, the sum is taken 
over all eigenstates, and (...) stands for averaging over 
the disorder. In contrast to the average density of states 



in conventional disordered systems the function P (e, y) 
at h 7^0 distinguishes between localized and extended 
states. For the localized states it is proportional to d (y) 
but is a non-trivial function for extended ones. 
It is convenient to rewrite Eq. (0) as 
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which allows to express the function P (e,y) in terms of a 
Gaussian integral over supervectors -0^ . This can be done 
using the fact that the ratio in the sum in Eq. (|j) can 
be composed of the elements of the matrix M^^ , where 



M. 
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In order to replace the integrals over all "0^ by integrals 
over supervector fields i/' (r) one should use vectors Uk (r) 
and Vk (r) and their conjugates Uk (r) and Vk (r) 
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where (t>k (r) and (f)k (r) are right and left eigenfunctions 
of the Hamiltonian H, Eq. g) (or Hl, Eq. (|)), the 
symbol * is complex conjugation. With the vectors Uk (r) 
and Vk (r) one can make Fourier expansion for any 2- 
component vector field. 

Using 8-component supervectors (r) with exactly the 
same structure as those in Rcfs. |1Q,P| one can rewrite 
Eq.(^ in terms of a functional integral with the La- 
grangian L 

C:^-i I i){v)ni!{r)dv (5) 

where the 8x8 matrix operator Ti has the form 

n = H' -e + i-ih + H"Ki + zyAiTa (6) 



InEq. 



H' = UH + H+) and H" = UH - H+) are 



Hermitian and anti-Hermitian parts of the Hamiltonian, 

respectively. In the continuum version H" = hV/m. 
The matrices A and t^ are the same as in Refs. [foP] 
and the matrix Ai anticommutcs with A being equal to 
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where 1 is the 4x4 unit matrix. A proper pre-exponential 
term is rather lengthy and is not written here. 

Further steps of derivation of the cr-model are stan- 
dard. One averages over the random potential and de- 
couples the effective interaction by integration over 8x8 
supermatrices Q. The integral over the eigenvalues of Q 



is calculated using the saddle-point approximation. The 
cr-model is finally obtained by expansion over VQ and h. 
As a result, the function P (e, y) takes the form 

2 /• 

P (e, y)^- lim ^ A[Q] exp {~F [Q]) DQ (8) 
where the free energy functional F [Q] is written as 

-4Q(7A + yAir3)]dr (9) 

In Eq. (||), Dq is the classical diffusion coefficient, [., .] 
is the commutator, STr is the supertrace, and v is the 
density of states of the system without disorder at h =0. 
The pre-exponential functional A [Q] equals 



A[Q] 



(Qi^(r) + gI^(r))(g^i(r') + Q^^(r')) 



- {Ql\ (r) + Q^i (r)) {Qll (r') + Qf^ {r'))]dvdv' (10) 

Numeration of matrix elements in Eq. (nO) is the same 
as in Refs. |l^,^; Eqs. (p|-|lC|) can be used in any dimen- 
sion. The supermatrix Q has the same symmetry as Q for 
the orthogonal ensemble in Refs. ||l^,^. The free energy 
functional F [Q] , Eq. (O) , has two additional terms with 
respect to the functional used for "conventional" disor- 
der problems. These terms contain the matrix Ai, which 
leads to new effective "external fields" in the free en- 
ergy. We see that the replacement of the physical vector- 
potential A by the imaginary quantity zh changes the 
symmetry of F [Q] . This refiects the fact that A and ih 
violate different physical symmetries. 

Of course, one can include in Eqs. (|l|,§) the vector- 
potential A, which would result in a standard term in 
Eq. (g) describing a magnetic field. If this field is strong 
enough one can use Eqs. (|-10) as before but in this case 
Q is a supermatrix with the structure corresponding to 
the unitary ensemble of Refs. [[lO|j9|]. This is a system 
with the broken time reversal invariance. 

Calculation of the function P (e, y) can be carried out 
using methods presented in Ref. p]. One can try to 
use renormalization group methods for study of 2D case 
(corresponding to 3D case for flux lines), transfer-matrix 
method for II? case (thick films with parallel line defects 
and magnetic field) or calculate definite integrals over Q 
for the QD (fiux lines in long cylinders). Leaving 11? and 
2D cases for future study let us consider now the QD case 
classified here as directed quantum chaos. 

In this limit one should integrate over Q assuming that 
this variable does not depend on coordinates. Then, the 
gradient terms in F [Q] in Eq. (0) can be omitted and 
integration over r, r' in Eq. (nol) easily performed. The 
QD form of F [Q] obtained in this way allows to make a 
very interesting conclusion even without starting explicit 



computation. This can be done comparing Eq. (||) with 
results of a recent work |l2] in which the supersymmetry 
technique was used to study density of complex eigenval- 
ues of "almost-Hermitian" random matrices X . These 
matrices were written in the form 



X = A + iaN-^/'^B 



(11) 



with N X N Hermitian statistically independent matri- 
ces A and B, and a number a of the order of unity. Due 
to the presence of the factor N~^^'^ the ensemble de- 
scribed by Eq. (O) differs from ensembles of matrices 
with arbitrary complex elements studied previously [[l3[ . 
Fyodorov et al |1^ calculated the joint probability of the 
real and imaginary parts of eigenvalues of the matrices 
X corresponding to the the function P (e, y), Eq. (g). 

Reducing in a standard way integration over the matri- 
ces A and B to integration over supervectors and then to 
supermatrices, the authors of Ref. |]1^ arrived in the limit 
A^ ^ cx) at the QD version of the free energy F [Q] , Eq. 
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The symmetry of Q corresponded to 



the supermatrices Q of Refs. ||10|,^ for the unitary ensem- 
ble, which is due to the hermiticity of the matrices A and 
B. Apparently, if they used real symmetric matrices A, 
antisymmetric matrices B, and imaginary a they would 
obtain Eq. (m with the supermatrices Q corresponding 
to the orthogonal ensemble (below these cases are called 
simply orthogonal and unitary). This demonstrates that 
the models of disorder, Eqs. (0,0) , in a limited volume are 
equivalent to the ensembles of weakly non-symmetric (or 
non-Hermitian if the time reversal symmetry is broken) 
random matrices. 

For explicit calculations, Fyodorov et al [Q used the 
parametrization of Ref. ||I^ . However, due to the pres- 
ence of the new terms with h and y in Eq. (0), the 
calculations with this parametrization are very difficult. 
So, the final result was obtained for the unitary ensem- 
ble only. As concerns the orthogonal case, computations 
with this parametrization do not seem to be possible at 
all. At the same time, study of the orthogonal ensemble 
can be very important because it describes the vortices 
in superconductors and, as will be seen later, the func- 
tion P (e, y) for the orthogonal ensemble at small h is 
qualitatively different from that for the unitary one. 

Fortunately, one can circumvent the difficulties using 
a new parametrization. Leaving details for another pub- 
lication I want to present here only a general structure. 
The supermatrix Q is written as 

Q - ZQoZ (12) 

with supermatrices Z satisfying the conditions ZZ — 1 
and [Z, Al] = 0. The central part Qo is chosen as 
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coscp 
-T3 sin( 



-T3 smt 
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The most lengthy part is calculation of the Jacobian but 
then the computation is of no difficulties and the density 
of complex eigenenergies P (e, y) for the unitary ensemble 
takes the form 



P e, y) = —IT- exp 



X 

1^ 



cosh xi exp {—at) dt 
(13) 



where x = 27ry/A, a? = 27riDo/i^/A, A is the mean level 
spacing. 

Eq. (|l3|) is exactly the result obtained in Rcf. ||l2[ . For 
the models of weak disorder, Eqs. (|lj||), P (e, y) depends 
only on the variable y because i^ is constant. For the 
random matrix model, Eq. (pi]), ly depends on e and 
obeys the Wigner semi-circle law. In the limit a 3> 1 one 
obtains 



Pi^,y) 



2a2A' 



< 2a^ 



(14) 



which means that the function P (e, y) is practically con- 
stant within the interval |a;| < 2a^. Beyond this interval 
the function exponentially decays. For the random ma- 
trix model this behavior corresponds to the well known 
"elliptic law" , Ref. [|l3[. For any a the function P (e, y) is 
smooth, which means, in particular, that the probability 
to find real eigenvalues is negligible. 

The situation drastically changes if the time-reversal 
symmetry is not broken. Although now the computa- 
tion of P (e, y) is somewhat more lengthy, one can finally 
express this function in a rather simple form 

P (e, y) = Pr+ Pc, Pr = lyS (y) / exp (-a'i') dt 

Jo 

Pc = ^$ 0-^) f exp {-ah^) sinh (|x| t) tdt 



(15) 



where $ (w) = ^ J^°° exp (— u^) du. 

We see that the function P (e, y) contains the part Pr 
proportional to 5 (y) . This means that a finite fraction 
of eigenstates has real eigenenergies for any finite a. The 
distribution of complex eigenenergies is described by the 
smooth function Pc. The fraction of the states with real 
eigenvalues vanishes only in the limit a —t oo. In this 
limit one comes for \x\ ^ a to Eq. (|l4|), which again 
results for the random matrix model in the elliptic law. 

The result about the finiteness of the fraction of eigen- 
states with real eigenenergies as well as Eq. (|l^) itself is 
new although indications of a peculiar behavior of the 
probability of real eigenvalues for asymmetric real matri- 
ces can be found in Ref. |Q. In Ref. ji) a mixture of 
states with real and complex eigenvalues was found nu- 
merically near the band center for the 2D model, Eq. (||) . 
Perhaps, the localization length in that region of param- 
eters exceeded the sample size. This would correspond to 



the 0-D situation rather than 2D. If it is so, the analytical 
result agrees with the numerical observation. 

Study of the directed disorder problems in higher di- 
mensions, although more difficult, is definitely of inter- 
est. Solution of the cr- model, Eqs. ( p^ , p^ ), in one di- 
mension seem to be possible using the transfer-matrix 
method. This case corresponds to a thick film in a par- 
allel magnetic field in the problem of the vortices in 
a superconductor. It would be interesting to obtain a 
localization-delocalization transition analytically. The 
one-dimensional cr-model is also known to describe the 
quantum kicked rotor problem ||l3|. May be, finite h cor- 
respond to a dissipation in this problem and one can find 
a transition here, too. 

In summary, quantum disorder problems with a di- 
rection (imaginary vector-potential) are considered. It is 
demonstrated that they can be mapped onto a new super- 
matrix cr-model. The QD version of the model describes a 
class of phenomena that can be called directed quantum 
chaos. Using the QD cr-model the equivalence of this 
class of systems to ensembles of weakly asymmetric or 
non-Hermitian random matrices is established. The joint 
probability of complex eigenvalues is computed explicitly 
and it is discovered that the fraction of real eigenenergies 
for time reversal invariant systems is always finite. 
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